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CLAIMS 



5 



1 . /vVmethod of evaluating assay data that Is arranged in a three 
dimensional array, the a^say being subject to systennatic and positional effects, 
the nnethod comprising: 

compensating the raw assay data for the systematic and 



10 positional effects; 



scoring the compensated data; and 

formatting the\(^ored data according to a determined format. 



2. The method of claim 1 wherein an assay is performed to 
1 5 generate a compendium of raw assay data that Is then compensated for 

systematic and positional effects. 

3. The method of claim 1 wherein the raw assay data is 
generated from a high throughput screening assay to identify a biologically active 

20 agent in a collection of test aCenfe, wherein a biologically active agent Is identified 
by identifying a test agent thargerferates a data point which statistically deviates 



from other data points in the 



ormatted scored data. 



4. The method of claim 1 wherein the raw assay data is 
25 generated from a high throughput screening assay to identify a pharmacologically 
active agent in a collection of test agents, 

wherein 

the high throughput screening assay is selected from the 
group consisting of enzyme activation, enzyme inhibition, llgand-receptor 
30 binding, ligand-receptor binding inhibition, cell cycle inhibition, cell cycle 
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activation, cell growth, cell division, cell activation, cell Inhibition, activation 
of production of and/or release of cellular factors, inhibition of production of 
and/or release of cellular factors, ion punnp, transport or channel activity, 
ion pump, transport or channel inhibition, activation of DNA synthesis, 
5 inhibition of DNA synthesis, activation of RNA synthesis, inhibition of RNA 

synthesis, activation of protein synthesis, inhibition of protein synthesis, 
nnetabollc activity, metabolic Inhibition, activation of apoptosis, and 
Inhibition of apoptosis, 

the collection of test agents is selected from the group 
10 consisting of: chemically synthesized molecules, natural products, cell 

extracts, nucleic acid molecules, cell culture media, proteins, isolated 
genetic material, fungal extracts and microbial fermentation broths; and 
recombinant products, viral particles, phage particles, proteins and peptide 
libraries; and 

1 5 the pharmacologically active agent Is identified by identifying 

a test agent that generates a data point which statistically deviates from 
other data points In the formatted scored data. 



5. The method of claim 1 comprising performing an assay such 
20 that a compendium of raw assay data is developed In a computer-readable form 

for purposes of allowing a computer-based algorithm to operate on such data. 

6. The method of claim 1 comprising compensating the raw 
assay data for systematic and/or positional effects by way of a computer-based 

25 algorithm. 

7. The method of claim 1 wherein the raw data has row-based 
aspects, column-based aspects, and non-additive Interaction-based aspects, the 
method comprising compensating the raw assay data for row-based positional 

30 effects. 
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8. The method of claim 1 wherein the raw data has row-based 
aspects, column-based aspects, and non-additive interaction-based aspects, the 
method comprising compensating the raw assay data for column-based positional 

5 effects. 

9. The method of claim 1 wherein the raw data has row-based 
aspects, column-based aspects, and non-additive interaction-based aspects, the 
method comprising compensating the raw assay data for longitudinal-based 

10 positional effects. 

1 0. A method of positionally correcting raw assay data from an 
assay comprising a plurality of longitudinally oriented plates p, each plate p 
having a plurality of wells organized into rows I and columns j, each well (i, j, p) 

15 having a raw value Xjjp associated therewith, the raw values x^p comprising the raw 
assay data, the method comprising deconstructing each raw value x^p of an 
associated well (i, j, p) Into: 

a plate effect value representing extraneous effects 
attributable to the plate p of the well (i, j, p); 
20 a row effect value representing extraneous effects 

attributable to the row I on the plate p of the well (i, j, p); 

a column effect value representing extraneous effects 
attributable to the column j on the plate p of the well (i, j, p); 

a non-additive. Interaction effect value representing 
25 extraneous positional effects attributable to consistent positional effects beyond 
the plate, row, and column effects previously determined for the (i, j, p) well on 
plate p; and 

a residual data value that is left over once ail the above 
extraneous effects are taken into account, 
30 the method further comprising employing the residual data value associated with 
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each well (i, j, p) to represent the well (I, j, p) as compared with all other wells (i, j, 
p) on the plate p. 

1 1 . The method of claim 1 0 wherein the raw assay data is 

5 generated from a high throughput screening assay to identify a biologically active 
agent in a collection of test agents, wherein the assay is subject to positional and 
systemic effects, the raw assay data is arranged in a three dimensional array, a 
biologically active agent is identified by identifying a test agent that generates a 
data point which statistically deviates from other data points in the formatted 
10 scored data. 

12. The method of claim 10 comprising employing the residual 
data value associated with each well (i, j, p) to represent the well (i, j, p) as 
compared with all other wells (i, j, p) on all of the plates p. 

15 

13. The method of claim 10 comprising deconstructing each raw 
value Xjjp of an associated well (i, j, p) into a plate effect value representing the 
overall median of all the raw values of the wells on plate p. 

20 14. The method of claim 13 comprising deconstmcting each raw 

value Xjjp of an associated well (i, j, p) into a row effect value representing the 
median of all the raw values of the wells for row i on plate p after taking into 
consideration the plate effect value. 

25 1 5. The method of claim 1 3 comprising deconstructing each raw 

value Xjjp of an associated well (i, j, p) into a column effect value representing the 
median of all the raw values of the wells for column j on plate p after taking into 
consideration the plate effect value. 
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16. The method of claim 10 comprising deconstructing each raw 
value Xjjp of an associated well (i, j, p) into a non-additive, interaction effect value 
representing an additional possible systematic measurement effect beyond the 
plate, row, and column effect values previously determined for the (i, j, p) well on 

5 plate p. 

17. A method of positionally correcting raw assay data from an 
assay comprising a plurality of longitudinally oriented plates p, each plate p 
having a plurality of wells organized into rows i and columns j, each well (i, j, p) 

10 having a raw value Xjjp associated therewith, the raw values Xj^p comprising the raw 
assay data, the method comprising: 

resistantly fitting the raw value x^p for each well (i, j) for each 
plate p to a row-column additive model: 

15 y^jp = Mp + R V + C'jp + e^p . 



where: 

y-p = the raw value for the well at row i and column j on plate p; 
Pp = an overall "average" for plate p; 
20 R'jp = a possible systematic measurement row offset for row i on 

plate p; 

C'jp = a possible systematic measurement column offset for column 
j on plate p; and 

ejjp = residual data without taking into account any non-additive 
25 interaction offset; 



each C'jp; 



longitudinally (plate-wise) non-linearly smoothing each RVp and 



substituting each un-smoothed R*jp and CVp value with a 
30 corresponding smoothed Rjp and C^p value and adjusting each residual value e^p 
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to an adjusted e'jjp: 

Yijp = Mp + Rip + + e'ijp ; 

5 longitudinally (plate-wise) non-linearly smoothing each e'^p to 

result in: 

Yijp = Mp + Rip + Cjp + smoothp(e'ijp) + r^p 

10 where each e'jjp is deconstructed into snnoothp(e'jjp), a possible systematic non- 
additive interaction offset for the (i, j) well on plate p, and r^p, residual data left 
over after taking into account any interaction offset; 

wherein each r^p represents a true relative value of the 
corresponding well (i, j, p) as compared to all other wells (i, j, p) on the plate p. 

15 

18. The method of claim 17 wherein the raw assay data is 
generated from a high throughput screening assay to identify a biologically active 
agent in a collection of test agents, wherein the assay is subject to positional and 
systemic effects, the raw assay data is arranged in a three dimensional array, a 

20 biologically active agent is identified by identifying a test agent that generates a 
data point which statistically deviates from other data points in the formatted 
scored data. 

19. The method of claim 1 7 comprising resistantly fitting the raw 
25 value Xjjp for each plate p to a row-column additive model according to a two way 

resistant median polish procedure. 

20. The method of claim 17 comprising longitudinally (plate-wise) 
non-linearly smoothing each R'ip and each C'jp according to one of a running 

30 median smoother and a lowess procedure. 
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21 . The method of claim 17 comprising longitudinally (plate-wise) 
non-linearly smoothing each e'^p according to one of a running median smoother 
and a lowess procedure. 

5 

22. The method of claim 17 further comprising normalizing each 
to result in a true relative value of the corresponding well (i, j, p) that can be 

compared to all other wells (i, j, p) on all plates p. 

10 23. The method of claim 22 comprising normalizing each r^jp by a 

standard deviation value derived from all the rijp's on the plate p to result in a 
score for the well (i, j, p) that can be compared across plates p: 

scorejjp = rjjp / (standard deviation value)p . 

15 

24. The method of claim 23 comprising normalizing each Tjjp by a 
median absolute deviation from median value. 

20 25. A computer having computer modules executing thereon for 

positionally correcting raw assay data from an assay comprising a plurality of 
longitudinally oriented plates p, each plate p having a plurality of wells organized 
into rows i and columns j, each well (i, j, p) having a raw value x^p associated 
therewith, the raw values x-^^ comprising the raw assay data, the modules 

25 comprising a first module deconstructing each raw value Xjjp of an associated well 
(i, j. p) Into: 

a plate effect value representing extraneous effects 
attributable to the plate p of the well (i, j, p); 

a row effect value representing extraneous effects 
30 attributable to the row i on the plate p of the well (i, j, p); 
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a column effect value representing extraneous effects 
attributable to the column j on the plate p of the well (i, j, p); 

a non-additive, interaction effect representing extraneous 
positional effects attributable to consistent positional effects beyond the plate, 
5 row, and column effects determined for the (i, j, p) well on plate p; and 

a residual data value that is left over once all the above 
extraneous effects are taken into account, the computer further comprising a 
second module employing the residual data value associated with each well (i, j, 
p) to represent the well (i, j, p) as compared with all other wells (i, j, p) on the plate 
10 p. 

26. The computer of claim 25 wherein the raw assay data is 
generated from a high throughput screening assay to identify a biologically active 
agent in a collection of test agents, wherein a biologically active agent is identified 

15 by Identifying a test agent that generates residual data value which statistically 
deviates from other residual data values generated. 

27. The computer of claim 25 wherein the second module 
employs the residual data value associated with each well (i, j, p) to represent the 

20 well (i, j, p) as compared with all other wells (i, j, p) on all of the plates p. 

28. The computer of claim 25 wherein the first module 
deconstructs each raw value Xjjp of an associated well (i, j, p) into a plate effect 
value representing the overall median of all the raw values of the wells on plate p. 

25 

29. The computer of claim 28 wherein the first module 
deconstructs each raw value x^p of an associated well (i, j, p) into a row effect 
value representing the median of all the raw values of the wells for row i on plate 
p after taking into consideration the plate effect value. 

30 
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30. The computer of claim 28 wherein the first module 
deconstructs each raw value Xjjp of an associated well (i, j, p) into a column effect 
value representing the median of all the raw values of the wells for column j on 
plate p after taking into consideration the plate effect value. 

31 . The computer of claim 25 wherein the first module 
deconstructs each raw value Xj^p of an associated well (i, j, p) into a non-additive 
interactive effect value representing possible systematic measurement effect for 
the (i, j) well on plate p beyond that attributable to the plate, row and column 
effect values. 

32. The computer of claim 25 further comprising an inputting 
module inputting the raw values Xjjp into a data structure in a memory of the 
computer, whereby the first module accesses the raw values Xjjp from the data 
structure. 

33. A computer having computer modules executing thereon for 
positionally correcting raw assay data from an assay comprising a plurality of 
longitudinally oriented plates p, each plate p having a plurality of wells organized 
into rows i and columns j, each well (i. j, p) having a raw value x^p associated 
therewith, the raw values x^p comprising the raw assay data, the modules 
comprising: 

a first module resistantly fitting the raw value x^p for each well 
(i, j) for each plate p to a row-column additive model: 

Yiip = Mp + R'ip + C' p + e^jp . 

the raw value for the well at row i and column j on plate p; 
an overall "average" for plate p; 



where: 

Yijp 
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R'ip = a possible systematic measurement row offset for row i on 
plate p; 

C'jp = a possible systematic measurement column offset for column 
j on plate p; and 

5 ejjp = residual data without taking into account any non-additive 

interaction offset; 



a second module longitudinally (plate-wise) non-linearly 
smoothing each R'jp and each C'jp; 
10 a third module substituting each un-smoothed RVp and C'jp 

value with a corresponding smoothed Rjp and Cjp value and adjusting each 
residual value ejjp to an adjusted e'^pi 



15 



20 



Yiip = Mp + Rip + Cjp + e'ijp ; and 

a fourth module longitudinally (plate-wise) non-linearly 
smoothing each e'jjp to result in: 



Yijp = Mp + Rip + Cjp + smoothp(e'ijp) + r. 



IIP 



where each e'jjp is deconstructed into smoothp(e'ijp), a possible systematic non- 
additive interaction offset for the (i, j) well on plate p and r^p residual data value left 
over after taking into account any non-additive interaction offset; 

wherein each r^p represents a true relative value of the 
25 corresponding well (i, j, p) as compared to all other wells (i, j, p) on the plate p. 



34. The computer of claim 33 wherein the raw assay data is 
generated from a high throughput screening assay to identify a biologically active 
agent in a collection of test agents, wherein a biologically active agent is identified 
30 by identifying a test agent that generates an which statistically deviates from r^p 




MERK-0004 / 20671 - 30 - PATENT 

generated. 

35. The computer of claim 33 wherein the first module resistantly 
fits the raw value x^p for each plate p to a row-column additive model according to 
5 a two way resistant median polish procedure. 



36. The computer of claim 33 wherein the second module 
longitudinally (plate-wise) non-linearly smoothes each RVp and each CVp according 
to one of a running median smoother and a lowess procedure. 



10 



J;=[ 37. The computer of claim 33 wherein the fourth module 



longitudinally (plate-wise) non-linearly smoothing each e'ijp according to one of a 
running median smoother and a lowess procedure. 

15 38. The computer of claim 33 further comprising a fifth module 

normalizing each rjjp to result in a true relative value of the corresponding well (i, j, 
p) that can be compared to all other wells (i. j, p) on all plates p. 

39. The computer of claim 38 wherein the fifth module normalizes 
20 each r^p by a standard deviation value derived from all the n^p's on the plate p to 
result in a score for the well (i, j, p) that can be compared across plates p: 

scorOiip = rfip / (standard deviation value)p . 



25 



40. The computer of claim 39 wherein the fifth module normalizes 
each rjjp by a median absolute deviation from median value. 



30 



41 . The computer of claim 33 further comprising an Inputting 
module inputting the raw values Xyp into a data structure In a memory of the 



MERK-0004 / 20671 



-31 - 



PATENT 



computer, whereby the first module accesses the raw values x^^ from the data 
structure. 

42. A computer-readable medium having computer-executable 
5 modules thereon for positlonally correcting raw assay data from an assay 

comprising a plurality of longitudinally oriented plates p, each plate p having a 
plurality of wells organized into rows 1 and columns j, each well (i, j, p) having a 
raw value Xjjp associated therewith, the raw values x^p comprising the raw assay 
data, the modules comprising a first module for deconstructing each raw value x^p 

10 of an associated well (i, j, p) Into: 

a plate effect value representing extraneous effects 
attributable to the plate p of the well (1, j, p); 

a row effect value representing extraneous effects 
attributable to the row 1 on the plate p of the well (1, j, p); 

15 a column effect value representing extraneous effects 

attributable to the column j on the plate p of the well (i, j, p); 

a non-additive, interaction effect representing extraneous 
positional effects attributable to systematic positional effects beyond the plate, 
row, and column effects previously determined for the (i, j, p) well on plate p; and 

20 a residual data value that is left over once all the above 

extraneous effects are taken into account, 

the computer further comprising a second module for employing the residual data 
value associated with each well (i, j, p) to represent the well (i, j, p) as compared 
with all other wells (i, j, p) on the plate p. 

25 

43. The computer-readable medium of claim 42 wherein the raw 
assay data is generated from a high throughput screening assay to identify a 
biologically active agent in a collection of test agents, wherein a biologically active 
agent is identified by identifying a test agent that generates a residual data value 

30 which statistically deviates from the other residual data value generated by the 
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assay. 

44. The computer-readable medium of claim 42 wherein the 
second module employs the residual data value associated with each well (i, j, p) 

5 to represent the well (i, j, p) as compared with all other wells (i, j, p) on all of the 
plates p. 

45. The computer-readable medium of claim 42 wherein the first 
module deconstmcts each raw value Xjjp of an associated well (i, j, p) into a plate 

1 0 effect value representing the overall median of all the raw values of the wells on 
plate p. 

46. The computer-readable medium of claim 45 wherein the first 
module deconstructs each raw value x^p of an associated well (i, j, p) into a row 

1 5 effect value representing the median of all the raw values of the wells for row i on 
plate p after taking into consideration the plate effect value. 

47. The computer-readable medium of claim 45 wherein the first 
module deconstructs each raw value Xjjp of an associated well (i, j, p) into a 

20 column effect value representing the median of all the raw values of the wells for 
column j on plate p after taking into consideration the plate effect value. 

48. The computer-readable medium of claim 42 wherein the first 
module deconstmcts each raw value x^p of an associated well (i, j, p) Into a non- 
25 additive interactive effect value representing a possible systematic measurement 

effect for the (i, j) well on plate p with respect to (i, j) wells on nearby plates p. 



49 The computer-readable medium of claim 42 further 
comprising an inputting module inputting the raw values Xjjp into a data structure in 
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a memory of a computer, whereby the first module accesses the raw values Xjjp 
from the data structure. 



50. A computer-readable medium having computer-executable 
5 modules thereon for positionally correcting raw assay data from an assay 
comprising a plurality of longitudinally oriented plates p, each plate p having a 
plurality of wells organized into rows i and columns j, each well (i, j, p) having a 
raw value x^p associated therewith, the raw values Xyp comprising the raw assay 
data, the modules comprising: 
1 0 a first module for resistantly fitting the raw value x^p for each 

well (i, j) for each plate p to a row-column additive model: 

yijp = Mp + Rip + C'p + eijp . 

15 where: 

y^p = the raw value for the well at row i and column j on plate p; 
Pp = an overall "average" for plate p; 

R'jp = a possible systematic measurement row offset for row i on 
plate p; 

20 C'jp = a possible systematic measurement column offset for column 

j on plate p; and 

eijp = residual data without taking into account any non-additive 
interaction offset; 

25 a second module for longitudinally (plate-wise) non-linearly 

smoothing each RVpand each C'^^\ 

a third module for substituting each un-smoothed R'jp and C'jp 
value with a corresponding smoothed Rjp and Cjp value and adjusting each 
residual value e-^p to an adjusted e'jjp: 

30 
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Yijp = Mp + Rip + + e ij, 



ijp » 



; and 



a fourth module for longitudinally (plate-wise) non-linearly 



smoothing each e'jjp to result in: 



5 



Yijp = ^Jp + Rip + C^p + smoothp(eijp) + Hj, 



ijp 



where each e'jjp is deconstructed into smoothp(e'ijp), a possible systematic non- 
additive interaction measurement offset for the (i, j) well on plate p, and rjjp, 
10 residual data left over after taking into account any interaction offset; 

wherein each r^p represents a true relative value of the 
corresponding well (i, j, p) as compared to all other wells (I, j, p) on the plate p. 



15 assay data is generated from a high throughput screening assay to identify a 

biologically active agent in a collection of test agents, wherein a biologically active 
agent is identified by identifying a test agent that generates an r^p which 
statistically deviates from the r^p generated by the other agents In the assay. 



module resistantly fits the raw value Xijp for each plate p to a row-column additive 
model according to a two way resistant median polish procedure. 

53. The computer-readable medium of claim 50 wherein the 
25 second module longitudinally (plate-wise) non-linearly smoothes each RVp and 

each C'jp according to one of a running median smoother and a lowess procedure. 



51 . The computer-readable medium of claim 50 wherein the raw 



20 



52. The computer-readable medium of claim 50 wherein the first 



30 



54. The computer-readable medium of claim 50 wherein the 
fourth module longitudinally (plate-wise) non-linearly smoothing each e'jjp 
according to one of a running median smoother and a lowess procedure. 
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55. The computer-readable medium of claim 50 further 
comprising a fifth module for normalizing each r^p to result in a true relative value 
of the corresponding well (i, j, p) that can be compared to all other wells (i, j, p) on 

5 all plates p. 

56. The computer-readable medium of claim 55 wherein the fifth 
module normalizes each r^^^ by a standard deviation value derived from all the n^p's 
on the plate p to result in a score for the well (i, j, p) that can be compared across 

10 plates p: 



module normalizes each r^p by a median absolute deviation from median value. 

58. The computer-readable medium of claim 50 further 
comprising an inputting module for inputting the raw values Xyp into a data 
20 structure in a memory of a computer, whereby the first module accesses the raw 
values Xiio from the data structure. 



scorejjp = rjjp / (standard deviation value)p . 



15 



57. The computer-readable medium of claim 56 wherein the fifth 



